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Abstract. The Karlsruhe Tritium Neutrino (KATRIN) experiment will measure the 
absolute mass scale of neutrinos with a sensitivity of m v = 200 meV/c 2 by high-precision 
spectroscopy close to the tritium /3-decay endpoint at 18.6 keV. Its Windowless Gaseous 
Tritium Source (WGTS) is a /3-decay source of high intensity (10 n /s) and stability, 
where high-purity molecular tritium at 30 K is circulated in a closed loop with a yearly 
throughput of 10 kg. To limit systematic effects the column density of the source has 
to be stabilised at the 10 -3 level. This requires extensive sensor instrumentation 
and dedicated control and monitoring systems for parameters such as the beam tube 
temperature, injection pressure, gas composition and others. Here we give an overview 
of these systems including a dedicated Laser-Raman system as well as several /3-decay 
activity monitors. We also report on results of the WGTS demonstrator and other 
large-scale test experiments giving proof-of-principle that all parameters relevant to 
the systematics can be controlled and monitored on the 10 -3 level or better. As 
a result of these works, the WGTS systematics can be controlled within stringent 
margins, enabling the KATRIN experiment to explore the neutrino mass scale with 
the design sensitivity. 
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1. Introduction 

With the discovery of massive neutrinos from //-oscillation experiments about a decade 
ago (see e.g. pQ and references therein), one of the most fundamental tasks for the 
next years will be the determination of the absolute mass scale of neutrinos. Neutrino 
masses can also be obtained from cosmological observations and neutrinoless double 
beta decay (e.g. [I]), but these methods strongly depend on the employed model. Model- 
independent results are provided only by kinematic /3-decay experiments [2J. 

The next-generation KArlsruhe TRItium Neutrino experiment (KATRIN) 
aims at improving the sensitivity in the neutrino mass measurement down to 
200 meV/c 2 (90%C.L.) [3], one order of magnitude lower than upper limits obtained 
by the Mainz [I] and Troitsk [5] experiments, the most sensitive direct neutrino mass 
experiments thus far. In this kind of experiment, information about the neutrino mass 
(i.e an upper limit or its value) was and will be derived with a precise measurement of 
the shape of the tritium /9-spectrum near its endpoint at E = (18571.8 ±1.2) eV |2J. 
The KATRIN experiment will be performed at the Karlsruhe Institute of Technology 
(KIT). Locating KATRIN at KIT allows the use of the unique expertise of the on-site 
Tritium Laboratory Karlsruhe (TLK), the only scientific laboratory equipped with a 
closed tritium loop [6] and licensed to handle the necessary amount of tritium required 
for the KATRIN experiment. 

The schematic layout of the KATRIN experiment is shown in figure [I] KATRIN 
intends to use a Windowless Gaseous Tritium Source (WGTS) with an activity of 
10 11 Bq. Such a source was employed by the LANL experiment [7] and developed further 
by the Troitsk experiment [8]. Superconducting solenoids will generate a homogeneous 
magnetic field of 3.6 T, which adiabatically guides the /3-decay electrons towards the 
tube ends. 

The energy measurement of the electrons is done by the spectrometer and detector 
section that is attached to the forward side of the WGTS via the transport section. The 
spectrometers are of MAC-E type [9j [TUJ. The detailed energy analysis is performed 
by the main spectrometer with an energy resolution of 0.93 eV. It uses an electrostatic 
retarding potential to transmit electrons with energies above the chosen retarding energy 
qU to the detector, while electrons with lower energies are rejected. By changing the 
retarding potential, KATRIN will measure an integrated tritium /3-spectrum, reaching 
its design sensitivity of 200 meV/c 2 (90% C.L.) within three full years of neutrino mass 
data taking, corresponding to about 5 to 6 years of operation. 
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Figure 1. Overview of the KATRIN Experiment. A more detailed description of the 
system components is found in the main text. 



The continuous precise knowledge or, even better, a high stability of the 
isotopic composition and the column density of the WGTS is very important for the 
determination of the neutrino mass. These parameters are associated with the main 
systematic uncertainties of KATRIN, namely activity fluctuations of the WGTS, energy 
loss corrections taking into account scattering of /3-electrons in the WGTS, and the final 
state spectrum [3J. For example, underestimating the energy loss will lead to a negative 
neutrino mass squared in the analysis, a problem that was encountered in analyses 
during the nineties of the last century [21 [TT]. Furthermore, the knowledge of these 
parameters is vital to combine data taken over the extended measurement period of a 
few years. 

Here we present the significant progress made in modelling the gas dynamics inside 
the WGTS as well as in the construction and commissioning of the actual components 
for controlling and monitoring the WGTS parameters. We focus on the methodologies 
and instrumentation needed to control and monitor the column density M and tritium 
purity €t as close to real time as possible, i.e. the time to obtain the information is 
shorter than or equal to the time interval (of the order of minutes) during which the 
retarding potential is constant before being shifted to the next value in the sequence. 
It should be noted that the different components are well past their design stage and in 
some cases are already fully operational. 

The paper is structured as follows: In section [2j we give an introduction on the 
WGTS and its important parameters. Our gas dynamics simulation is presented in 
section |3j as well as the resulting requirements for the different control and monitoring 
methods for the WGTS. Their status is described in section |4} The paper finishes with 
a presentation of the impact for KATRIN in section [5| 

2. The Windowless Gaseous Tritium Source (WGTS) 



Key parameters of the WGTS to achieve high statistical sensitivity in the neutrino 
mass search are the column density Af, defined as the number of molecules within a flux 
tube volume of unit cross section (further discussed in section 2.1.1), and the isotopic 
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composition (T 2 , DT, HT, D 2 , HD, H 2 ), as the count rate S of the source scales like: 
S = C ■ e T ■ AT. (1) 



Here refers to the tritium purity further discussed in section 2.1.2; the proportionality 
constant C encompasses various experimental properties (detector efficiency, acceptance, 
etc.). 

It should be noted, that equation is not necessarily valid over the full TV- 
range since the observed activity of a source saturates for large column densities. The 
reason for this is that electrons escape the source only with reduced energy, due to the 
increasing scattering probabilities, which is enhanced by the increasing cross section 
towards lower energies; thus these electrons cannot overcome the retarding potential 
anymore. This can lead to non-linearities in the relation between the measured source 
activity and column density. However, simulations with full particle tracking of the 
scattered electrons in the magnetic field down to a threshold of 100 eV show that the 
non-linearity of equation ([I]) is negligible in an interval ±1% around a central value 
Af = 5.0 x 10 17 molecules/cm 2 . 

Both parameters, Af and ex, also play a significant role for the systematic 
uncertainty in deriving m^from the measurement as TV strongly affects the energy 
losses of the electrons inside the source by inelastic scattering. Non-T 2 impurities lead 
to mandatory corrections of the /3-electron spectrum, because of their different recoil 
energies and different molecular final state distributions. 

In the discussion of the requirements for Af and €t monitoring, we will distinguish 
between trueness, precision and accuracy (e.g. [T2]). The trueness A tme is defined as 
the difference between the mean value of a measured observable and the true value, 
for example A tme (A0 = (A/" mea s) — Aftme- The precision A prec corresponds to the 
square root of the statistical variance of several measured values around their mean 
value, for instance A prcc (Af) = a/ ( (Af mcas - ( A/" meas ) ) 2 ) • The accuracy A acc is defined 
as the difference between a single measured value and the true physical value, e.g. 

A acc (A/") Af mcas Af true- 

In section 2.1 the entities Af and e-r are defined in detail together with their targets 
in trueness and precision. Those are derived from the KATRIN goal, namely that no 
single systematic effect gives rise to an uncertainty of Am 2 < 7.5 • 10~ 3 eV 2 /c 4 to the 
neutrino mass analysis [3]. In section 2.2 the technical realization of the WGTS is 
presented with the focus on the required control and monitoring systems to access Af 
and €t experimentally. 



2.1. Key Parameters of the WGTS 

From a physicist's point of view, the WGTS is a column of tritium gas inside an 
open cylinder with a diameter of D = 90 mm and a length of L = 10 m inside 
a homogeneous magnetic field of B = 3.6 T generated by superconducting solenoid 
magnets (see figure [2]). The gas consists of different species with mole fractions Cj. It is 
injected near the middle with an adjustable pressure pi n of about 10~ 3 mbar, and pumped 
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Figure 2. Principle of the Windowless Gaseous Tritium Source (WGTS). The density 
profile inside the beam tube is kept constant by continuously injecting tritium gas in 
the middle and pumping out at both ends. 



out at both ends with constant pumping speed with an outlet pressure p out ~ 0.05 • p m . 
Apart from the geometry and this pressure gradient, the flow for a gas with viscosity 
r] is determined by the beam tube temperature, which has a certain profile T B t(<P, z) 
with a mean value in the 30 K regime. Here, z denotes the cylinder axis with z = 
being the middle in the source, the azimuthal angle. 

These parameters define the distribution function fi(r,v) describing the number of 
molecules of species i, which can be found in the phase space volume [r, r + d 3 r], [v, v + 
d 3 v] . From calculations of /$ (r , v) , the influence of the above parameters on the column 
density M can be investigated. Details on this are given in section |3.1 



The full computation of the distribution function fi(r,v) is also required for 
investigations of plasma and ion physics effects (e.g. calculations of space charge 
effects |14j ) or for investigating the interactions of low-energy secondary or shake-off 
electrons due to the short path length of such electrons. Also, the velocity distribution of 
the mother molecules is determined by the aforementioned fi(r, v). Thus one may apply 
thermal Doppler-broadening corrections to the /3-electrons associated with the thermal 
diffusion and position-dependent gas flow speeds of the decaying tritium molecules (see 
section 



3.1). Such corrections are essential for the exact description of the /3-electron 



spectrum. 



2.1.1. Column density The column density M is defined as 

+L/2 

N(Pm,Pout,T B T((fi, z)) = J ni(p,(fi,z,p in ,p ou t,T BT ((j),z))dz (2) 

* -L/2 

with rii(p, 0, z,p in ,p out ,T BT ((p, z)) being the integral of /j(r, v) over the velocity phase 
space in cylindrical polar coordinates (p, <f), z) for a given set of experimental boundary 
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conditions. In fact M is of course a function of (p, <fi), but due to the homogeneity in 



(p, (j>) (see section 3.2) this dependence is omitted in the notation. 



At first glance, the knowledge of the true value A/true seems of no particular interest 
as the signal strength and therefore J\f (see equation ([I])) is an unrestricted parameter 
in the fitting algorithms of the neutrino mass analysis as long as activity fluctuations 
during one scan of the /3-spectrum are monitored and taken into account. However, 
the neutrino mass analysis depends on the accurate description of inelastic scattering of 
the signal electrons by the gas molecules inside the source (cr iriel (18.6 keV) = 3.40(7) • 
10 -18 cm 2 ) [13J. Under nominal conditions, 41.3% of all electrons escape the source 
without inelastic interaction, 29.3% are scattered once and 16.7% are scattered twice. 
Hereby, the reference column density of the WGTS of M = 5.0 x 10 17 molecules/cm 2 and 
the nominal 5-field configuration of KATRIN, which accepts /3-electrons with starting 
angles with respect to the S-field of up to 51° for transmission through the KATRIN 
beam line, is assumed. In principle, these scattering probabilities and associated energy 
losses will be measured directly in a dedicated calibration measurement [3]. However, 
such a measurement is time-consuming and thus cannot be applied in parallel to the 
neutrino mass search. Hence, monitoring the amount of energy losses due to scattering is 
required. This is achieved by continuously measuring M with an experimental precision 
A prec (A0. 

According to neutrino mass sensitivity simulations, the necessary precision 
A wec (N)/N to keep the corresponding systematic error below Am 2 < 7.5 • 
10~ 3 eV 2 /c 4 [3] depends on details of the neutrino mass analysis, in particular the width 
of the analyzed energy window below the endpoint. For reference, we use the most 
demanding value from [3J, namely A piec (J\f) / J\f < 0.1%. 

As can be seen from equation ([2]), changes in J\f are due to shifts or fluctuations of 
the experimental conditions of inlet and outlet pressures Ap in , Ap out and temperature 
changes AT(r, 0). The latter can either be localized or affect the whole beam tube. 



2.1.2. Isotopologue composition The tritium purity €t is defined as the ratio of the 
number of tritium atoms to the total sum of atoms in the WGTS. A high tritium purity 
is primarily necessary to maximize the signal rate 5* according to equation ([!]). It is 
planned to maintain €t > 0.95 over the complete measurement period. This will be 
achieved using molecular tritium (T 2 ) as the main gas constituent with a high mole 
fraction c(T 2 ) > 0.9, a small admixture of DT (c(DT) < 0.1) and only trace amounts of 
HT, D 2 and H 2 . The tritiated hydrogen isotopologues T 2 , DT and HT differ not only in 
their mass and therefore in their respective recoil energies, but also in the rotational and 
vibrational final state distribution of their daughter molecules following tritium decay. 
Both effects are taken into account in the modelling of the WGTS. The corresponding 
final state distribution for T 2 for the first four initial rotational states and DT/HT for 
the initial states J = and 1 have been calculated in [TS] . 

The systematic uncertainty for m 2 that is introduced by the theoretical description 
of the different final states is estimated to be 6 ■ 10~ 3 eV 2 /c 4 [3J. In the neutrino mass 
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analysis, the final state distributions of the individual isotopologues will be weighted 
with the measured mole fraction of the individual isotopologue. The accuracy of the gas 
composition measurement therefore has to be converted into a systematic uncertainty 
A acc (gas) for m 2 which has to be combined with the uncertainty of the final state 
description, i.e. (A acc (m 2 )) 2 = (6 • 10~ 3 eV 2 /c 4 ) 2 + (A acc (gas)) 2 . The uncertainty 
A acc (gas) should be below a certain limit (1 • 10~ 3 eV 2 /c 4 ) in order to avoid a significant 
inflation of the combined uncertainty even in case of an improvement of the theoretical 
description of the final state distributions. Simulations [16J show that an accuracy of 
A acc (eT) ~ 1% is required, which can also be derived from general considerations [2]. 
This is a rather moderate requirement due to the envisaged high T 2 mole fraction 
(e T > 0.95). 

Regardless of this, there is compelling reason to measure with at least a 
magnitude better precision. Namely, it opens up a complementary approach for 
monitoring A (A/") which is independent of hydrodynamical calculations of /j. According 
to equation ([!]), monitoring M is equivalent to monitoring both and the activity of 
the WGTS. For this reason the WGTS instrumentation is extended by different activity 
detectors and a Laser Raman spectroscopy system (LARA) to measure with a relative 
precision of A prec (eT)/eT < 0.1%. 



2.2. The WGTS cryostat 

Achieving the proposed tritium purity and activity in the WGTS constitutes a rather 
challenging task. Specifically, the column density inside the source beam tube has to be 
kept constant to a level of at least 0.1 %. Accordingly, a variety of parameters have to 
be stabilized; these are: 

- the injection rate, 

- the pumping speed in forward and rear direction, 

- the gas composition, 

- and the beam tube temperature. 

Fulfilling these requirements results in a complex cryostat-system, the so-called WGTS 
cryostat [T7|. Its cryogenic system consists of 13 fluid circuits operated with 6 cryogenic 
fluids. The cryostat requires a measurement and control system with more than 500 
sensors. Besides the actual WGTS, the cryostat incorporates the DPS1-F and DPS1-R 
(the first stages of the differential pumping sections in the front and rear direction, 
respectively). Also included are the first seven superconducting magnets guiding the 
electrons adiabatically from the source beam tube through the transport section to the 
spectrometers. 

The DPS1-R will be connected to the Calibration and Monitoring System (CMS) 
via the so-called rear-wall. The rear-wall separates the WGTS beam tube from the CMS 
and defines the electrical potential of the tritium plasma inside the WGTS. Furthermore, 



it may be used to measure the /3-electron current directly (see section 4.4.1 ), provided the 
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Figure 3. The cooling concept of the WGTS beam tube is based on a thermo-siphon 
with two 16 mm evaporator tubes brazed onto the beam tube, which are partially filled 
with neon at its boiling point (30 K). Convection in the thermo-siphon will be controlled 
by heater wires inside the evaporator tubes. The vapour flows to the evaporator ends, 
and back through a common return line to the condenser. The thermo-siphon is cooled 
with liquid helium. 



current measurement can be shielded from low energetic secondary electrons and ions. In 
addition, it acts as an X-Ray converter which allows activity measurements of the source 
via Beta-Induced X-Ray Spectroscopy (BIXS, see section 4.4.2). Furthermore, the CMS 
houses an e-gun for calibration measurements (see section 4.5). In the forward direction, 
the DPS1-F is connected to the DPS2-F, a further differential pumping section. Its role 
is to reduce the tritium flow by at least another four orders of magnitude [18J and to 
host an FT-ICR (Fourier Transform Ion cyclotron resonance) mass spectrometer for ion 
physics [19] . 

The basic concept of keeping the injection rate into the WGTS constant is based 
on injecting tritium gas from a pressure controlled buffer vessel over a 5 m long capillary 
with constant conductance and pumping the gas at both ends of the beam tube. The 
details of the complex tritium gas handling system, the so-called Inner Loop, are 
discussed in section 14. 1L 

The Inner Loop is supplied with high purity tritium by the Outer Loop System, 
which is covered by the existing infrastructure of TLK [6J. The isotopic purity in the 
Inner Loop System is monitored in-line with a Laser Raman Spectroscopy setup (LARA) 
discussed in detail in section 14.31 

The set-point for the WGTS temperature is about 30 K. This constitutes a near 
optimum: at higher temperatures, the gas flow rate and the throughput increases. Thus 
a larger T2 inventory would be required to reach the same column density; also the 
Doppler broadening effect increases. At lower temperatures, clustering and condensation 
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of the hydrogen isotopologues would be induced which would be equally undesirable. 

The concept of providing the cooling for the 10 m long source tube is summarized 
in figure [3| As shown there, it is based on a two phase neon thermo-siphon. Two 
copper tubes are brazed over the whole length on the outside of the source beam 
tube and thus determine the temperature on both sides. The neon gas flows to a 
condenser in a closed cycle; there it is liquefied again, being cooled with gaseous 
helium at ~25 K. The measured fluctuations of the helium cryogenic circuit used 
to cool the condenser are of the order of 300 mK per hour [2T]. These fluctuations 
exceed the temperature stability requirements derived in the next section. Therefore, 
these fluctuations need to be damped by the coupling between the helium circuit and 
the two phase neon system, therefore the condenser has a large heat capacity. Since 
such a cooling concept has never been applied on a large component like the WGTS, 
a full scale test experiment (the WGTS demonstrator) has been performed in 2011. 
This demonstrator consists of original components of the final WGTS assembly, in 
particular the beam tube and its instrumentation. The main difference to the final 
WGTS cryostat is that all superconducting magnet components are replaced by a cold 
mass to simulate their thermal properties during the demonstrator tests. Furthermore, 
the WGTS demonstrator was operated without external potential while the WGTS in 
its final configuration can be operated on a variable voltage of up to 1000 V to define 
the potential difference between WGTS and main spectrometer. Selected results of the 



demonstrator experiment are briefly summarized in section 4.2 and will be discussed in 
more detail in a future publication [20] . 

The aim of the demonstrator tests is to understand the thermal behaviour of the 
system and to verify that the only non-negligible influx of heat into the beam tube 
originates from the pump ports at both ends of the WGTS tube, as predicted and 
shown in |21j . This heat influx leads to an azimuthal temperature profile of the beam 
tube that can be measured by temperature sensors distributed along the tube and on 
several azimuthal positions. The measured temperature distribution can then be used 
in the gas dynamics modelling for the WGTS. 

3. Modelling of the WGTS gas dynamics 



The first part of this chapter, section |33j focuses on the current state of the gas dynamics 
simulation developed specifically for the WGTS to compute the density profile and 
column density, as well as the velocity profile of molecules in the source. From these, one 
can determine the influence of experimental parameters on the column density and thus 
ultimately on the neutrino mass sensitivity. In addition, one can deduce the specification 
for the monitoring of the crucial WGTS parameters presented in section |3.2| 
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3.1. The gas dynamical model 

In order to determine the influence of the thermodynamical properties of the WGTS on 
the column density J\f as defined by equation (|2]), a detailed gas dynamical simulation 
of the WGTS has been developed, which allows one to compute the density distribution 
n(r) and consequently Af. The velocity distribution f(p,<f),z,v) is also of interest, 
since the resulting electron energy is Doppler-shifted due to the motion of the tritium 
molecules. 

One can characterize the flow regime in the WGTS beam tube with the so-called 
rarefaction parameter 



S = w ith v m = J — —. (3) 

r] v m V m 

Here R is the tube radius, p the pressure, rj the viscosity and v m is the most probable 

speed for molecules with mass m and temperature T. Dependent on the flow regime 

one can distinguish three cases: 

- The hydrodynamic regime with 5 ^> 1, where the equations of continuum mechanics 
are valid. 

- The free molecular flow regime (5 <C 1), where intermolecular collisions can be 
neglected and the molecules move independently of each other. 

- The transition regime, where the full Boltzmann equation needs to be solved, using 
a detailed model of intermolecular collisions. Since for this the intrinsic nature of 
the intermolecular interactions is important, this is the most difficult scenario. 

The rarefaction parameter varies from 5 ~ 20 at the injection in the WGTS beam tube 
to 5 < 10 _1 at the beginning of the differential pumping chambers, with a significant 
part of the beam tube exhibiting values within the intermediate regime. One can use 
the phenomenological intermediate conductance formula of Knudsen [23] to compute 
the gas distribution approximately inside the WGTS, delivering an easily computable 
formula for the pressure distribution in [24]. To be more precise and to be able to 
include detailed boundary conditions, in principle the full Boltzmann equation needs to 
be solved. 

The length of the beam tube L is much larger than its radius R. Thus, to a good 
approximation, the radial dimensions can be neglected (n(f)= n(z)) as it will be shown 



in subsection 3.2 The problem is reduced to only a single (the longitudinal) spatial 
dimension, which simplifies the analytical description and tremendously reduces the 
necessary CPU time required for the calculation. 

For this one-demensional case, the method outlined in [25] can be applied: It relies 
on the fact that the mass flow rate 



M 



7TR 



■3 



-G P (6)-f + G T (5) 



(4) 



dz T(z) dz 

driven by small gradients of pressure p and temperature T, has to remain constant 
along the tube. The Poiseuille coefficient Gp and the thermal creep coefficient Gt, 
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Figure 4. (a): Calculated relative density profile n(z)/n(0) along the WGTS 
beam based on an expected source temperature distribution, an injection pressure 
of Pin = 3.368 • 10~ 3 mbar at z = and exit pressures of p out — 0.05 p m at 
z = 5 m. (b): The expected relative deviations from the average column density 
(M(x, y) — A/"avc) /Wave for a quarter of the WGTS cross section and an assumed 
azimuthal temperature variation. 



which depend only on the rarefaction parameter 5, are available in [22] ■ The coefficients 
were obtained by linearising the Boltzmann equation around a Maxwell-Boltzmann 
distribution, employing the S-model j2S] for intermolecular collisions and using the small 
pressure and temperature gradients as perturbation parameters. 

With these coefficients, and incorporating the boundary conditions T(z) and p in (cf. 
sections 4.1 and 4.2), provided by the continuous monitoring systems, equation Q can 
be solved numerically. As a result one can determine the pressure p = p(z), and density 
n(z) = p(z)/kBT(z) profiles along the tube. The resulting density profile obtained for 
typical values of these boundary conditions is shown in the left part of figure |4| 

However, the one-dimensional model cannot fully describe the gas flow in the 
more complex geometries of the injection and pump port region. Furthermore, the 
injection and exit pressures cannot be measured directly, but have to be inferred from 
pressure sensors in the Inner Loop. Therefore, the development of more detailed three- 
dimensional simulations of these regions is currently in progress. As a first step, a radial 
density distribution caused by an azimuthal temperature variation was computed for 
different positions along the beam tube. Weighting these cross sections with their tritium 
content allows the computation of the expected radial column density fluctuations shown 
in the right part of figure |4| 

It should be stressed that the gas dynamics simulation will not be used to 
quantitatively determine the absolute value of M. Although it is, in principle, possible 
to compute the column density for a given set of boundary conditions T(z),p in and 
p out , the simulation has only a trueness of about 5%. The trueness is limited by two 
facts: Firstly, important fluid dynamical properties like the viscosity of tritium and 
the accomodation coefficient between tritium and the beam tube walls are not known 
accurately enough. Secondly, intrinsic uncertainties of the simulation, in particular the 
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Table 1. Proportionality constants between different experimental parameters and 
the column density [3J. 



parameter 


a x 




temperature 


C(t = 


-1.2 


injection pressure 


a Pin = 


= 1.1 


exit pressure 




= 0.03 



modelling of the intermolecular collisions and using the linearized Boltzmann equation to 
compute the coefficients Gp and Gt in equation Q, also cause systematic uncertainties 
within the simulation. Solving the exact Boltzmann equation numerically has only been 
achieved in very few and simple cases; the complex WGTS case would require excessive 
computational effort and is thus not pursued. Instead, the gas dynamics simulation is 
used to determine the influence of a certain (experimental) parameter on TV. Once this 
influence is known, the task of monitoring M is translated into continuously measuring 
this parameter with the required precision. In particular the existing one-dimensional 
simulation is sufficiently fast to be suitable as a nearly real-time analysis tool for the 
column density based on the actual real-time monitoring of the beam tube temperature 
and injection pressure. 



3.2. Results of the gas dynamics simulation 

Of particular interest for the monitoring requirements of the WGTS is the influence of an 
experimental parameter X = T ,p in ,p out on the column density defined in equation 
which governs the energy loss of /3-electrons in the source: 
A/V AX 

This influence has been determined with the gas dynamics simulation in [3] and is 
summarized here. The coefficients ax given in table [T] allow for determining the stability 
requirements for the corresponding physical quantity. 

With the current version of the gas dynamics simulation, it has become possible to 
include the influence of a temperature profile on n(f) and M too, as determined by the 



instrumentation on the WGTS beam tube (see section 4.2). Now, it is also possible to 
compute the full distribution function f(z,v) instead of only the density distribution. 
It turns out that the velocity distribution of the tritium molecules is well described 
analytically using a so-called local Maxwell-Boltzmann distribution 

f{p,<f>,z,v) = exp . (6) 

W nV m) \ V m J 

The bulk velocity U z (p, <fr, z) reflects the fact that the gas flows from the injection in the 
middle of the beam tube to both ends. U z (p, 0, z) is of the order of 10 m/s for most of the 
beam tube length, but increases to about 50 m/s at the beam tube exits. This is much 
smaller than the most probable thermal speed v m = \j2kpT '/nix ~ 288 m/s, where nix 
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is the molecular mass of the respective isotopologue. According to the Doppler effect, 
the motion of the decaying molecule changes the /^-electron energy as 



For example, the additional energy change AE for an 18.6 keV /3-electron emitted by a 
tritium molecule with velocity v m corresponds to AE = cos9re • 130 meV, where 9re is 
the angle between momenta of the mother molecule and the /3-electron. Since this is of 
the order of the KATRIN neutrino mass sensitivity, it needs to be taken into account 
in the /3-spectrum calculation. 

Several additional investigations have been carried out to study specific effects that 
potentially influence the column density: 

• A detailed simulation of the injection region based on [27J has been performed to 
optimize the design of the injection chamber; a design with many small holes was 
shown to be advantageous to avoid turbulence. The final design of the injection 
chamber features 415 holes with a diameter of 2 mm. 

• The thermal velocity of the different hydrogen isotopologues differs slightly due to 



their mass difference yv^ = i/mT 2 /mDTfm ~ l-l^m 2 j- This results in separation 
phenomena and consequently a dependence of the gas composition on the position. 
This was studied in dedicated simulations based on |28j . In the standard operational 
mode of the WGTS, these effects are negligible (the relative change is below 10 -4 ). 
However, in a special operational mode of the WGTS, the so-called Krypton Mode, 
this effect becomes significant. In the Krypton Mode, a small fraction of 83 Kr is 
added to the gas mixture for calibration purposes. This implies that the WGTS has 
to be operated at 120 K since krypton freezes at lower temperatures. In this case, 
the larger temperature in combination with the larger mass of krypton increases 
de-mixture effects significantly, but this is not relevant for the work described here. 

• The influence of an azimuthal temperature variation was investigated and can be 
calculated for inhomogeneities up to a few K by solving the linearized Boltzmann 
equation using the discrete velocity method outlined e.g. in [22} [291 ED]- However, 
as shown in figure [4], this effect is negligible for the expected radial temperature 
gradients since the corresponding column density fluctuations are below 10~ 4 . 

Overall, according to these simulations, the influence of the corresponding effects on 
n(r) and therefore A^ is of the order of AM /M < 10~ 4 ; this is safely below the critical 
value for an unaccounted shift of AJ\f crit /Af = 0.1% [3]. Consequently, monitoring 
T(z),p in ,p ont is sufficient. 

The effect on the shape of the tritium /3-spectrum by tritium decays occurring in 
the pump ports or in the differential pumping section (DPS1-F, DPS1-R) has not yet 
been investigated in detail. In |3], it has been noted that effects from decays outside 
the central beam tube should in principle not be significant due to their small fraction 
(Pout < 0.05p in ). To properly account decays from these regions, which are partially 
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Table 2. Reference values and requirements for the stability of WGTS parameters, 
partially from J3J. 



parameter 


reference value 


req. relative stability req. trueness 


temperature 


T = 


(27- 30)K 


10~ 3 


0.5 K 


injection pressure 


Pin = 


= 3.368 • lCr 3 mbar 


10~ 3 




exit pressure 


Pout 


* (0.03 - 0.05)p in 


0.03 




purity 


e T > 


0.95 


10~ 3 


0.1 


activity 


A = 


10 n Bq 


10~ 3 





operated at 80 K and which have inhomogeneous magnetic fields, a simulation with full 
3-dimensional geometry has to be done in the future. 

The resulting stability requirements for the monitoring systems are summarized in 
table [2j In addition to the precision requirements obtained with the values from table [TJ 
the absolute value of the beam tube temperature is also of interest; it influences the 
integrated /3-spectrum through the Doppler effect. However, the impact of the Doppler 
effect is relatively small due to the low source temperature. In recent investigations, we 
have shown that the resulting accuracy requirement is only about 500 mK. 

For completeness, we briefly address the requirements for gas composition and 



activity monitoring. As outlined in section 2.1.1, monitoring these two quantities 
constitutes a complementary approach to the monitoring of Pi n ,Pout and T(z), since 
column density fluctuations can also be obtained by measuring both the gas composition 
and the source activity using equation Q. Furthermore, as a consequence of the 
different final state distributions of T 2 , DT and HT, now also the trueness (and not 
only the precision) of the gas composition measurement is of relevance. Fortunately, the 
resulting accuracy requirement is much less stringent than for the precision; in recent 
studies, we have shown that in case of the envisaged high tritium purity > 0.95, a 
1% accuracy for ex can be reached with a trueness of the LARA system of only about 



4. Control and monitoring 

In table [2] above, the actual requirements for the control and monitoring of the crucial 
source parameters are listed. These requirements can only be met if the gas inlet and 
outlet rates, the beam tube temperature and the gas composition are stable. In the 
Technical Design Report of KATRIN [3j only rough concepts were sketched. Now that 
research and development has been performed on all relevant control and monitoring 
devices, it has been demonstrated that these requirements can be fully met. 
The different control and monitoring systems summarized in figure [5] are: 

- The Inner Loop system provides stable gas injection and the stable gas 
composition together with the infrastructure of the Tritium Laboratory Karlsruhe 



(TLK) (see section 4.1). 
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- A sensor system of PT500 temperature sensors provides the monitoring of 
the spatial and temporal temperature distribution T BT (<p,z) on the WGTS beam 
tube (see section 4.2). 

- The Laser Raman system (LARA) is connected to the Inner Loop and enables 
the in-line monitoring of gas composition (see section 4.3). 

- Three activity detectors are foreseen to measure the activity of the WGTS. The 
overall WGTS activity is proportional to the product of column density (Af) and 
isotopic purity (e^) (see section 4.4). 

- An electron gun embedded into the Calibration and Monitoring System (CMS) 
of the KATRIN experiment is used to measure the scattering losses of the electrons 



in the WGTS (see section 4.5). 



It should be noted that most of the systems now have been set up, and we have proven 
their suitability for the various control and monitoring tasks. The few remaining systems 
are already beyond the prototype phase and the results thus far suggest that they could 
be implemented into KATRIN in the near future. Below we introduce the various control 
and monitoring systems, and discuss their current status. 



4-1. The inner tritium loop 

From the point of view of tritium circulation, achieving a stable column density requires 
(i) the stable injection of tritium in the middle of the WGTS beam tube, and (ii) stable 
pumping on both ends. 

While conceptually this sounds very simple, in reality it requires a rather complex 
assembly. Overall, the actual working principle very much resembles Ohm's law (for 
gases): Q = Ap ■ c. A constant pressure gradient Ap and a constant conductance c 
(which will be achieved by a constant temperature and length of the connection pipe) 
will result in a constant gas flow Q into the WGTS. 
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Figure 5. Overview of the analytic systems for control and monitoring of 
the gaseous tritium source. The detailed description of the systems is found in the 
main text. Note that both spectrometers are scaled to 50% of their size. 
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Figure 6. (a) Simplified flow scheme and (b) picture of complete setup of the tritium 
circulation loop. 

To fulfil the first requirement, a highly precise "flow controller" had to be built at 
the Tritium Laboratory Karlsruhe; this is the so-called Inner Loop (figure [6j (b)). It 
consists of 45 sensors for pressure, temperature and ionization, 70 automatic and manual 
valves, 2 flow controllers, 1 flow meter, 9 vacuum pumps, 6 buffer vessels, 2 palladium 
membrane filters, 70 custom made pipes, ~ 20 metal bellows and a total of > 400 
metallic seals. 

The second requirement, i.e. a stable pumping rate, will be achieved using turbo 
molecular pumps (Oerlikon Leybold MAG W 2800) on both ends of the WGTS. Due 
to a constant rotation frequency (on a 10~ 4 level), a constant pumping speed and 
consequently a constant exit pressure will be maintained. 

After setting up the flow controller, leak tests of every single component and an 
integral leak test were performed. The measured integral leakage rate of 7- 10~ 9 mbard/s 
conforms with the regulations for primary systems of the Tritium Laboratory. 

Since the WGTS is not yet available for test measurements, the connection had to 
be bypassed by a tube with similar conductance as the later combination of WGTS and 
connection line. 

An example for the measured pressure trend during pressure-stabilized circulation 
of deuterium is shown in figure [7j Here, the KATRIN limits are indicated by the dotted 
lines; clearly, the pressure stabilization falls well within these limits. Note that the 
particular run shown here may be seen as the worst-case. The apparent dip during the 
early times of the run did not occur in any other of the runs; it could well have been 
the result of a variation in the room temperature. From the trace in figure [7j (a), one 
extracts an average pressure (p\) = 15.024 ± 0.002 mbar inside the pressure-controlled 
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Figure 7. (a) Results of an individual pressure-control run. The pressure in the 
WGTS buffer vessel is shown for several hours of deuterium circulation and a pressure 
set-point of 15 mbar. The dashed line represents the ltr requirement for KATRIN 
operation, (b) Measured pressure fluctuation as a function of the mean pressure during 
nitrogen circulation, as measured for different set-points. 



buffer vessel, for a set-point of 15 mbar. 

The corresponding gas flow is stable within 120.4±0.5 seem, with a temperature 
fluctuation of the buffer vessel of ±0.125 K at 318 K. Fourier analysis of the pressure 
trends shows no indication for periodic short term (seconds/hours) or long term 
(days/weeks) fluctuations. 



For calibration (see e.g. section 4.5) and verification reasons, the pressure 



stabilization is designed to operate between pi ^ 1 mbar and pi ~ 20 mbar. A 
series of pressure stability tests at different pressures has been made with N 2 filling. 
The results are shown in figure [7] (b). The observed fluctuations are well within 
the KATRIN requirement for all pressure set-points. In particular, for the standard 
KATRIN conditions with a pressure set-point of p\ ~ 15 mbar and only negligible 
temperature fluctuations on the connection line to the WGTS, the Inner Loop performs 
nearly an order of magnitude better than the 0.1% limit required by KATRIN, for both 
deuterium and N 2 . 

To examine the effect of a temperature variation of the WGTS connection line, the 
WGTS bypass is equipped with an electrical heater to simulate temperature changes 
resulting in a different conductance. The temperature T bypass was varied in four steps 
over an interval of « 20 K within a few hours. This exceeds the expected fluctuations 
under KATRIN conditions by three orders of magnitude. Such extreme variations of 
up to 20 K do not affect the pressure to such an extent that the KATRIN 0.1% limit 
would be surpassed. For a variation ATb yP ass of ~ 2 K (still ~ 100 times the KATRIN 
specification), no visible influence on the pressure stabilization has been observed. 
Therefore, it can be concluded from these test measurements that the Inner Loop fulfils 
and even exceeds the requirements for KATRIN. 
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4-2. Beam tube temperature 

In both the WGTS demonstrator setup and the final WGTS assembly described in 



section 2.2 the monitoring of the beam tube temperature is performed by 24 metallic 
resistance temperature sensors (Pt500) and 24 vapour pressure sensors. The Pt500 
sensors are used to monitor the temperature continuously Their resistance is measured 
by digital multimeters with a resolution of AR = 3 ■ 10~ 5 Q. A common calibration 
curve, measured for the original sensors of the WGTS, is used to convert the resistance 
to actual temperatures with a theoretical temperature resolution of AT = 3 • 10 -5 K. 
Several contributions to the uncertainty have been evaluated in |3T] ; the results suggest 
a combined calibration uncertainty of 125 mK at T = 30 K; the main contribution 
seems due to the magnetic field dependence of the Pt500 temperature measurements. 
This uncertainty is insufficient to reliably monitor the beam tube temperature with 
requirements of 30 mK on stability and homogeneity. Therefore, each Pt500 sensor 
is calibrated in situ by an adjacent vapour pressure sensor. A vapour pressure sensor 
is a small volume filled with LNe. The temperature is determined by measuring the 
saturation pressure of neon with a transducer outside of the cryostat, connected to the 
sensor by a small capillary. This measurement is insensitive to magnetic fields and 
has a measurement uncertainty of only 4 mK. However, continuous measurement are 
not possible since the neon filling of the sensors can neither be guaranteed nor verified 
during the expected KATRIN run periods of up to 60 days. Instead, the very accurate 
vapour pressure sensors can be used to calibrate the Pt500 sensors before and after 
every KATRIN tritium run to reduce the uncertainty of each individual sensor to that 
of the LNe pressure sensors, i.e. 4 mK. This is sufficient to monitor the beam tube 
temperature. 



As mentioned in section |2.2[ this approach was tested in 2011 with the WGTS 
demonstrator test setup. Data on the stability of the primary LHe supply and 
thermosiphon (see figure [3]) have already been analysed. 

In particular, the shielding of the two-phase neon system from the mass flow 
and temperature fluctuations of the primary helium circuit, which are of the order 
of 300 mK [21], works as expected. This is evident from the discrete Fourier transforms 
of both the primary helium circuit and beam tube temperature shown in figure [8] While 
the helium circuit exhibits significant fluctuations and characteristic frequencies, only 
uncritical low-frequency fluctuations on the time scale of days are transferred from 
the helium circuit through the condenser to the two-phase neon system; all critical 
frequencies with / > 0.3 mHz are reduced by at least two orders of magnitude. This 
indicates that the temperature stability is well within the design specifications for 
KATRIN. However, the full analysis of the data taken at the WGTS demonstrator 
is still ongoing and will be reported in a future publication 
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Figure 8. (a) Discrete Fourier transform of the temperature of the original 
cryogenic helium used to cool the neon system for a 2 day measurement period, (b) 
Corresponding Fourier transform of the signal of one of the temperature sensors on 
the beam tube for the same period. 



4-3. Gas composition monitoring using Raman spectroscopy 



As stated in section 2.1.2 , the isotopic purity of the inlet gas of the WGTS has to be 
monitored continuously with a precision < 10~ 3 . The integration time for an individual 
measurement should be as short as possible to be compatible with the measurement 
intervals of KATRIN. Typical methods for quantitative analysis of tritium, e.g. gas 
chromatography [32] and liquid scintillation counting [33], are not suitable for non- 
stop operation at the KATRIN experiment because of the duration of the measurement 
procedure (typically 3> 250 s), the need of sample extraction and the continuous waste 
production. In contrast, Laser Raman spectroscopy is a non-contact and in-line method 
which is based on the inelastic Raman scattering [31] of light from gas molecules 
and continuous monitoring in measurement intervals of < 250 s is feasible. Raman 
spectroscopy systems are commercially available and widely used for gas analysis. But 
when combining the requirements of operation with tritium and < 10 -3 statistical 
uncertainty, off-the-shelf commercial products are not suitable [35] . 

Although, in principle, no other gas species than hydrogen isotopologues are 
expected to be injected into the source tube, Raman spectroscopy is also sensitive to 
any other molecular species like nitrogen and tritiated methane species (CT 4 , CDT 3 , 
etc.) [36]; the former is thought to originate from incomplete evacuation of the loop 
prior to the filling with tritium, while the latter can be attributed to exchange reactions 
with carbon liberated from the steel vessel and tubes. 

Over the past five years, a specialized Laser Raman spectroscopy unit (LARA) has 
evolved and was thoroughly tested [371 EH]- It consists of a measurement cell (LARA 
cell), which is located between the buffer vessels of the Inner Loop (see figure [6J and a 
laser setup for gas analysis (see figure |9j (a)). The gas pressure inside the LARA cell 
is typically in the range of 150-200 mbar during operation of the Inner Loop. A DPSS 
laser, operating at 532 nm, is used as excitation source. The Raman scattered light is 
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Figure 9. (a) Simplified schematic drawing of the LARA setup (for details see main 
text). The polarization angle of the laser beam is controlled by a half- wave plate (A/2). 
The photo diode and power meter are used for monitoring of the laser stability. The 
LARA cell is located inside the glove box of the Inner Loop while the optical setup 
is on the outside to prevent contamination, (b) Representative Raman spectrum of a 
gas mixture with p to t = 217 mbar and er > 90%. The Raman branches are labeled 
according to [31]. T 2 dominates the mixture but traces of DT and HT are visible. 



collected under an angle of 90° by optical lenses and guided by an optical fibre to a 
spectrometer (Princeton Instruments, HTS) and CCD detector (Princeton Instruments, 
PIXIS:2KB) for spectral analysis. The Raman shifts and intensities of the vibrational 
Qx branches of the hydrogen isotopologues are used for qualitative and quantitative 
analysis, respectively (see figure [9j (b)). A custom-written analysis package has been 
developed, that incorporates the automated correction of high intensity CCD pixels due 
to cosmic rays, background subtraction and peak fitting [39J. 

Although the present system is still based on the originally proposed concept [3], 
recently three major improvements have been implemented which aided in maximizing 
the signal-to-noise ratio: 

- A temperature controlled high-stability laser (Laser Quantum, Finesse 5 W, 
532 nm), and thermally stable mounts for the optical components are used to 
reduce fluctuations of laser power and beam alignment that would impose intensity 
fluctuations in the Raman spectra. 

- The laser beam is passed through the LARA cell twice to nearly double the effective 
laser power in the cell with respect to the standard, single-pass setup. 

- Groups of CCD pixels are simultaneously read-out ("on-chip binning") to minimize 
noise from the analogue-to-digital conversion in the CCD |39j . 

By 2011, the LARA system was tested in the tritium test circulation loop, LOOPINO. 
Here it had already reached A prec (eT 2 )/eT 2 — 0.1% for T 2 monitoring within 250 s 
acquisition time under conditions similar to KATRIN operation (~200 mbar absolute 
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pressure, tritium purity €t ~ 0.95) before the aforementioned improvements were 
implemented [36]. At this time, a level of detection of 0.012 mbar partial-pressure 
equivalent was reached [35J. After implementation of the improvements, a first test with 
a tritium mixture (p(Tg) ~ 7.4 mbar, p(HT) ~ 1.6 mbar, p(DT) rj 0.7 mbar) was 
performed. For this, a statistical uncertainty of A prec (eT 2 )/eT 2 = 0.3% was achieved 
for T2 monitoring within 250 s acquisition time. Extrapolating this performance to 
the KATRIN operating conditions, the LARA system can reach A prec (eT 2 )/eT 2 = 0.1% 



within 60 s acquisition time. 

In contrast to the other monitoring tools of the WGTS, not only the precision 
but also the trueness of LARA measurements is of interest, due to the influence of the 
different final state distribution of the individual isotopologues on the /3-spectrum (see 
section 2.1.2). 

To address the problem of trueness in the observed Raman signal amplitudes, 
a calibration method with non-trititated reference gas mixtures has been developed 
to determine the response functions for H2, HD and D2, which relate the Raman 
line intensities to the molar fraction of the isotopologues in the gas mixture. The 
response functions, which are different for the individual isotopologues, include quantum 
mechanical transition matrix elements, the z/ 4 wave number dependence of Raman 
scattering and the spectral sensitivity of the light collection system of LARA. In this 
study, an apparent variation between the individual response functions by up to 10% 
has been observed. 

In an alternative approach, "synthetic" calibration spectra are generated against 
which the measured spectra are compared. For this, the transition matrix elements of 
the hydrogen isotopologues are calculated using ah initio formalisms [10] • Unfortunately, 
the relevant matrix elements were not yet experimentally verified for all hydrogen 
isotopologues; hence, any slight errors therein would potentially propagate into the 
LARA measurement of the gas composition. Therefore, a precise experimental check 
of the transition matrix elements of all hydrogen isotopologues has been performed and 
results will be reported shortly [4"Tj . 

Another calibration issue in determining the trueness in the Raman peak 
amplitudes is the spectral efficiency of the light collection system; for this a complete 
characterization of the LARA setup is required, which will rely on a well-defined 
calibration light source (presently under development). The aim is to further improve 
the accuracy of the LARA calibration (currently ~ 10%). However, should the 
understanding of the final states improve, and thus the uncertainty estimate reduce, 
then this would immediately impact on the LARA accuracy requirements. Most likely 
they would become more stringent; however, the current improved implementation of 
the LARA setup would certainly be capable to meet such a demand. 

As a final note, we would like to mention that during operation with concentrated 
tritium gas mixtures, damage of the tritium facing anti-reflection coatings of the LARA 
cell windows was observed [SB]. A test experiment is currently ongoing to search for 
tritium resistant coatings and to investigate the potential loss of optical transparency of 
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the windows under /3-radiation. For final testing, a long-term measurement campaign 
of the LARA setup with a circulating tritium gas mixture (ex > 0.95) is planned for 
2012. 



4-4- Activity measurements 

The activity monitoring systems presented in this section have two important goals: (i) 
They provide information about fluctuations of the WGTS activity itself with a time- 
scale of minutes; these would cause systematic effects in the neutrino mass analysis if not 
taken into account, (ii) Together with the measured tritium purity ey, equation ([!]) can 
be used (in the form Af = Si/ (Ci ■ e^)) to monitor fluctuations of the column density J\f 
with 0.1 % precision. Here, the index % reflects the fact that the proportionality constant 
Ci and thus the signal rate Si is different for each activity monitoring system. The 
combined monitoring of both activity and gas composition thus allows one to disentangle 
whether observed activity fluctuations are due to changes of the gas composition or 
changes in the column densitity A/"(which have a different effect on the KATRIN count 
rate). 

Various concepts to measure the source activity have been studied. In all methods, 
a lower threshold for the /3-electron detection is necessary to separate the signal from 
the overwhelming amount of low-energy secondary and shake-off electrons. Furthermore, 
the detection method must not disturb the electron transport from the source to the 
main spectrometer. 

In the following, we present three experimental concepts and devices which can 
be integrated into the KATRIN beam line. Two of these are located in the Control 
and Monitoring System (CMS), with either directly measuring the current in the 



rear wall incurred by the /3-electrons (see section 4.4.1), or with Beta-Induced X-ray 



spectroscopy(see section 4.4.2). The third device, the Forward Beam Monitor Detector 
(see section 4.4.3), is located in the Cryogenic Pumping Section (CPS). Prototype tests 
demonstrated that with all three concepts the detection of electrons with a suitably low 
energy threshold is feasible. 

The column density will also be measured periodically (every few hours) by an 



electron gun; as discussed in section |4.5[ this is an important calibration measurement 
by itself. It also allows to probe whether the constants Ci of the activity monitoring 
systems are truly time independent (within the required precision of 0.1%). 



4-4- 1 ■ The rear wall as a Faraday cup To use the rear- wall as a Faraday cup to monitor 
the WGTS source activity, the WGTS /3-electron emission is treated as a DC current, 
which is of the order —10 nA. This is accompanied by low-energy plasma currents, which 
must be screened. For this a multi-layer thin-film coating is applied to the rear wall. The 
top layer is a gold film of about 50 nm thickness; while relatively thin, this gold layer is 
sufficient to intercept all low-energy plasma species. It is electrically isolated from the 
metallic substrate by a 150 nm Si02 insulator of 10 6 — 10 7 Q. This dual-layer thin-film 



Monitoring of the operating parameters of the KATRIN WGTS 



23 





-4 _ o Au 



o 



o 



-5 







2 



_ 4 6 8 
Electron energy /keV 



10 



12 



Figure 10. (a) Cross-section through a multi-layer rear-wall sample, (b) Measurement 
of the currents from the aluminium and gold layers in a different sample, as a function 
of the electron beam energy. This measurement displays the ability of a multilayer 
rear wall to separate low-energy from high-energy electrons. 



structure is sufficient to meet the aforementioned threshold conditions: electrons above 
~5 keV reach the substrate, while lower-energy electrons are stopped in the gold layer. 
A high-precision ammeter can then measure the substrate current to 0.1% precision in 
0.1 s. The radiation dose to the Si02 layer is of the order of 20 MGy/y. In fact at such 
doses performance degradations for Si0 2 capacitors are not expected. Such exposures 
will be tested in advance also to check the stability of the space charge accumulation. 

A mix of atomic-layer deposition (ALD) and plasma-enhanced chemical vapour 
deposition (PECVD) on an Al test substrate, is used to produce the Au/Si0 2 thin-film 
structures described above. For aluminium test substrates we achieved intact insulator 
areas as large as 30 cm 2 . These test samples were exposed to an electron beam of a 
few nA, and the current impinging onto the Au-layer and through the Al-substrate were 
measured simultaneously, as a function of electron energy. One observes that low-energy 
beams stop in the gold, high-energy beams penetrate down in the substrate, with the 
crossover energy (equal currents in both layers) at ~ 6 keV. The cross section of the 
layer structure for a selected sample is shown in figure [TOj together with preliminary 
current measurements at different electron energies. The displayed data verify that this 
multilayer design fulfils the anticipated, basic functionality. Further systematic studies 
are under way (i) to measure the uniformity of the gold properties; (ii)to expose the 
insulator to high radiation doses and (iii) to produce a sample of full size (200 cm 2 ) on 
the desired beryllium substrate. 

If the on-going R&D is successful and if it is possible in such a current 
measurement to discriminate the arriving primary /3-electrons from the low-energy 
secondary electrons, then this method will allow to monitor column density changes 
over the whole source cross section within measurement times of just a few seconds. 

4-4-2. X-ray detection for measuring the source activity Most of the electrons from 
tritium /3-decay in the WGTS (« 10 11 per second) are magnetically guided to the rear 
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Figure 11. Schematic of the TriReX setup. 



wall (including the electrons originally flying towards the spectrometer section, as most 
of them are reflected by the retarding potential and eventually hit the rear wall). The 
absorption of /3-decay electrons in the gold coating of the rear wall then generates X- 



rays (see section 2.2), which can be detected using Beta-Induced X-Ray Spectroscopy 
(BIXS). 

Provided that the rear wall substrate is sufficiently transparent for X-rays produced 
by bremsstrahlung or fluorescence, a detection system for activity monitoring of the 
WGTS can be placed behind the rear wall. 

A proof of principle experiment (Tritium Rear wall experiment - TriReX) was built 
at TLK to demonstrate the feasibility of monitoring the activity of gaseous tritium under 
conditions similar to the WGTS. 



A simplified setup of TriReX is shown in figure [TT] Conceptually, it is divided into 
three separate volumes. The first (right) recipient volume contains the (H,D,T)-gas 
mixture. An axial magnetic field up to 120 mT provided by a water-cooled air coil 
guides the /3-decay electrons to the gold coated beryllium windows at both ends. The 
right beryllium window (of 50 mm diameter and 180 /im thickness) serves as a reference 
surface for later adsorption measurements, while the left beryllium window (with 39 mm 
diameter and 200 /im thickness) constitutes the equivalent of the KATRIN rear wall. In 
addition, the window prevents contamination of the silicon drift detector (SDD, KETEK 
AXAS-M), located in the second (central) volume, with tritium. The SDD (active area: 
80 mm 2 ) is oriented towards the beryllium window. The energy resolution of the system 
is less than 160 eV at 5.9 keV. The third (left) volume contains the detector electronics, 
including pre-amplification and temperature control. 

A spectrum taken with the TriReX setup is shown in figure 12, (a). The continuous 
bremsstrahlung spectrum is superimposed by X-ray fluorescence lines. The X-ray line 
spectrum is dominated by a peak around 2 keV. This peak originates from the Zr- 
collimator which is factory mounted to the SDD. Additional peaks can be associated 
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Figure 12. (a) Bremsstrahlung spectrum and characteristic X-rays, measured with 
the TriRex setup, together with the calculated bremsstrahlung spectrum for a tritium 
partial pressure of (0.8 ± 0.1) mbar and a total activity of (1.4 ± 0.1) • 10 11 Bq. The 
total intensity of the calculated spectrum was fitted to the measurement, (b) Corrected 
integral count rate vs. tritium partial pressure (error bars are included) taken with 
the TriReX setup. The corrected integral count rate is the integral count rate minus 
the memory effect signal N = (35.2 ± 0.3) cps after the evacuation of the system. The 
starting pressure of 0.098 mbar equates to (8.3 ±0.8) • 10 9 Bq. The system was pumped 
down in successive steps to a final pressure of 0.001 mbar. 



with Cr, Fe, Mn and Mo which are present as trace constituents in the stainless steel 
vacuum vessel. Two further peaks at higher energies are caused by the gold coatings of 
the beryllium windows. The characteristic X-rays provide the unique possibility of an in- 
situ energy calibration of the system. The cut-off below 2 keV is caused by absorption of 
the X-rays within the beryllium windows; the energy threshold of the detector had been 
set to < 150 eV. The shape of the measured bremsstrahlung spectrum agrees well with 
calculations using the analytical model from jl2] and the mass attenuation coefficients 
taken from |43j. 

With the observed count rate of N = 14486 cps, the necessary statistical precision 
of A prec = 0.1 % is reached in t < 70s measurement time within TriReX. For pressures 
lower than ^10 mbar the integral count rate depends in principle linearly on the activity 
in the recipient [H]. Deviations from linearity as seen in figure 12 (b) are most likely 
caused by outgassing effects. The TriReX setup could not be baked out because of 
design limitations. Therefore outgassing from the walls during the pump-down causes 
an activity signal which is lower than expected because of the decreasing tritium purity. 
Also a tritium memory effect has been observed because of tritium adsorption on the 
walls. This memory effect causes an increasing background signal after each tritium 
measurement run. 

The achievable count rate in the Control and Monitoring Section (CMS) for the 
WGTS setup depends on (i) the final geometry of the rear wall and the X-ray detector; 
(ii) the thickness and the layer composition of the rear wall; and (iii) the number of 
implemented X-ray detectors. However, the measured spectra here agree well with the 
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calculated and simulated expectations. Therefore, a reliable extrapolation to the WGTS 
case is viable. 

For the activity monitoring of the WGTS, the current design of the CMS includes 
two SDDs, with 100 mm 2 active area each. The thickness and the layer composition 
of the rear wall is still under investigation with R&D efforts to coat a ~ 400 /im 
beryllium window with 120 nm of mono-crystalline gold. The expected count rate for 
this particular material and coating configuration would be ~ 6.6 kcps, allowing the 
activity monitoring with 0.1 % precision within a sampling time of ~ 150 s. 

Further investigations regarding the long term stability of such a system, 
backgrounds by tritium adsorption on gold and stainless steel and an absolute calibration 
of the system at the ~ 1% level will take place in the near future. 



4-4-3. Forward Beam Monitor Detector Silicon radiation detectors, which directly 
detect the /3-electrons in the flux tube, can also be used to monitor the source activity. 
Such a detector should not shadow any parts of the electron flux tube, which is used for 
the neutrino mass analysis, neither in the rear nor the forward directions. Therefore, 
care has to be taken where to position the detector. Because of these constraints, 
this component will be implemented in the so-called Forward Beam Monitor Detector 
(FBMD) configuration, which probes only the outermost rim of the flux tube. A 



picture of the setup is shown in figure 13 It is assumed that this localized activity 



measurement is nevertheless representative for the whole beam cross section. That this 
assumption holds may be verified during periodic calibration runs, since the detector 
concept also includes the option to move the detector element across the whole flux 
tube cross section. The FBMD will be mounted between the last two superconducting 
solenoids of the cryogenic pumping section (CPS), i.e. at the end of the transport section, 
where the Tritium partial pressure is reduced by 14 orders of magnitudes compared to 
the WGTS. This minimizes any background effects and contamination of the sensitive 
detector surface by tritium, and vice versa minimizes any impact from the detector on 
the Tritium gas. 

1.6 T and can 
e with a position 



B 



The magnetic field B at this position is axially symmetric with 
therefore be used to study the spatial homogeneity of the source profi 
resolution of the manipulator of 50 /im. For its standard monitoring role, the detector 
will be 61-67 mm away from the beam axis. 

The concept was tested by a proof-of-principle experiment using an off-the-shelf 
PIN diode (manufactured by Micron semiconductor). It comprises a circular active 
area of 1.3 mm 2 , a thickness of 300 /im and a dead layer of approximately 100 nm. The 
diode chip is mounted with epoxy on a ceramic substrate which also carries the charge- 
sensitive pre-amplifier as a hybrid module. The signal path is split into a spectroscopic 
AC mode and a DC mode which can be used in parallel. Both methods can be utilised 
to measure the actual activity of the WGTS. 

In the AC mode, the energy of each impinging electron is measured individually by 
digitalization of the pre-amplifier signal and subsequent pulse height analysis. 
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Figure 13. Setup of the Forward Beam Monitor Detector, (a) Manipulator with 
detector module in indicated flux tube (same scale), (b) PIN diode with an active area 
of 1.3 mm 2 . 

In the DC mode, the continuous flow of decay electrons results in an increased 
DC offset of the integrating pre-amplifier, which scales with the activity. Of course, 
information on the energy of the /5-electrons is lost in the DC mode, but sampling times 
are significantly reduced with respect to the AC mode; thus the activity can be measured 
with 0.1% precision within just a few seconds. 

The reason for this is that in DC mode no "pile- up " effects are encountered, which 
are prevalent at higher count rates in AC mode. As a consequence, a larger detector 
area could be chosen for the DC mode detection, capturing a larger electron flux. The 
opposite would hold for the AC mode: count rates can be limited by making the detector 
area smaller. Currently, the sensitive detector area is a compromise between the optimal 
0.1mm 2 for AC readout mode and 5 mm 2 which would be the optimum for the DC 
readout mode. For the final design, it is proposed to use a custom designed detector 
chip with two active pixels, their respective sizes optimized for AC and DC readout 
mode. 

For both modes, temperature stabilization of the detector is important. In the 
spectroscopic (AC) mode, the energy resolution and thus the stability of detection 
threshold depends on the detector temperature. In the DC mode, the diode leakage 
current which is superimposed on the signal produced by the impinging electrons, rises 
exponentially with the detector temperature. Therefore, a precise temperature control 
and stabilization of the detector element at temperatures of about — 30°C is necessary. 
The cooling concept is based on a highly flexible copper braid with a cross section of 
10 mm 2 which thermally connects the movable detector holder to a copper heat sink. 
The heat sink is cooled by a cold gas supply which is based on liquid nitrogen in a 
dewar vessel evaporated with an electrical heater. With this approach, temperature 
stabilization to better than ±0.1 K has been achieved. 

The detector as well as its electronics and temperature control was mounted 
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on a two dimensional UHV manipulator; the complete assembly was tested with an 
electron gun. The beam scanning capability of the manipulator was successfully 
shown. In the (spectroscopic) AC mode the monitoring of the electron flux proved 
to be reliable up to rates of about 50 kHz when pile-up becomes the limiting factor. 
Within the linear response regime, the electron energy resolution was measured to be 
cx(FWHM) = 1120±50eV@17.5keV; this agrees with calculated noise sources like Fano 
and electronic noise. 

In the (integral) DC mode, count rates of up to 212 kHz have been measured with a 
systematic uncertainty of 1.9%. It should be noted that during long-term measurements 
drifts of the DC signal were observed; it was suspected that electrons hitting the corners 
of the detector surface create this effect. A detector optimized for this measurement 
mode, incorporating an additional shielding, should reduce or avoid these adverse effects. 

It was possible to show that the FBMD system is capable to provide activity 
monitoring to a reasonable precision, however for the final KATRIN implementation 
some modifications are required. The dimensions at the KATRIN experiment dictate 
that a larger UHV manipulator, but of the same type and the same cooling concept, 
has to be used. The detector, especially its sensitive area, needs to be designed to 
accommodate the high /3-intensity of 1.6 • 10 6 s _1 mm -2 . Due to the low energy of the 
electrons, which is insufficient to create lattice defects, radiation damage is not an issue 
[35] . But the detector sensitive area needs to be limited to reach an impinging rate of 
e.g. 50 kHz in the AC mode, meaning that within an integration time of 30 s an activity 
uncertainty of 0.1% can be achieved. On the other hand, the DC detector can measure 
with a high rate like 5 MHz to reach a precision of 0.1% within 1 s integration time, 
assuming that the systematic error remains constant. 

4-5. Other measurement methods for the column density 

In this paper, the focus has been on tools for in-line monitoring and characterization 
of the WGTS, on time scales as close as possible to real-time. For completeness, 
we briefly address a complementary, but off-line, method to infer A/", namely via the 
analysis of scattered electrons. The method is based on an approach which involves 
both the high-resolution main spectrometer and the focal plane detector. Consequently, 
the precision of the method is linked to the performance and precision of those two 
main KATRIN components. It has the potential to directly measure the energy losses 
caused by scattering, with a nominal accuracy of the probability for n-fold scattering 
P n of A acc (P n )/P„ < 0.1%, and it may be seen as one of the most crucial calibration 
measurements for the KATRIN experiment. 

Unfortunately, this monitoring method cannot be performed simultaneously with 
the neutrino mass measurements. Therefore, it will be used only in periodic tests rather 
than for continuous monitoring. Despite not being a near real-time method it is a vital 
cross-calibration check for Af and therefore it is worthwhile to briefly describe it here. 

An electron emitter (electron gun, e-gun) in the CMS based on the photo-electric 
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effect produces quasi-mono-energetic electrons with stable intensity and with a well- 
defined angle 9 (the angle between electron momentum and 5-field) and adjustable 
kinetic energy Eq [16]. The parameter pair (Eq,6) can be varied over the full range of 
tritium- /3-electron values encountered in the WGTS, with a maximum angular spread of 
oq < 4°. The maximum kinetic energy is technically limited to 25 keV with a proposed 
maximum energy spread of <te < 0.2 eV in the low rate mode (R< 10 4 cps). 

The measurement principle relies on the correlation between the electron scattering 
probability and column density M : changes in the WGTS column density can be invoked 
by the changes of probabilities for no, single or multiple scattering events. The scattering 
probabilities are solely dependent on the total cross section, the path length through 
the source and the column density. From the Poisson distribution, one finds that the 
probabilities for no-loss and single scattering of an electron entering the WGTS from 
the CMS amount to: 

F o (0) = e"^ an d P X {Q) = ■ (8) 

cos 

Due to the high energy resolution of the main spectrometer, the energy loss 
spectrum &N/&E\ oss can be measured for a given condition of the electron gun (E ,9). 
Then, deconvolving the differential cross section da/dEi oss from the measured spectrum 
yields P n [3]. The required differential energy loss function is either taken from models 
or derived from a special calibration measurement. According to the formulae (|8j), the 
column density M can be monitored with the measured precision of A(Pi/P ). However, 
the method is not suitable as near-time monitoring since this would lead to a significant 
loss of measuring time for the neutrino mass campaign. Therefore, it is useful only 
for periodic examinations instead of continuous monitoring. Nevertheless, it can be 
regarded as complementary method to the other AJ\f measurements described above. 

5. Conclusions 

In order to achieve the design sensitivity in the measurement of the neutrino 
mass, all main-task components in the KATRIN experiment need to operate with 
extreme stability over extended periods of time, and their behaviour needs to be 
well characterized. In particular, this is true for the WGTS which provides the 
essential tritium /3-electrons from whose energy spectrum analysis the neutrino mass 
will be extracted. Here, we have presented the status of the control and monitoring 
system that is designed to ensure that the stability of the key WGTS parameters, 
namely the column density A" and tritium purity ey can be monitored and analysed 
as close to real-time as possible. For this task, theoretical modelling and experimental 
methodologies were developed; both were successfully combined confirming that the 
challenging requirements for the WGTS operation can be met. 

Gas dynamics simulations were used to link the requirements for the precision 
and accuracy of A" and to experimentally observable quantities like pressures and 
temperatures. Progress in this modelling includes (i) the incorporation of temperature 
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inhomogeneities in the gas dynamics simulation, and (ii) a detailed investigation of the 
position-dependent velocity distribution in the WGTS. 

The first strategy to monitor M relies on monitoring the experimental parameters 
influencing A/", namely injection and outlet pressure as well as the beam tube 
temperature, with a precision determined by these simulations. The pressure stability is 
realized by the Inner Loop system. The Inner Loop is operational and has reached and 
even exceeded the necessary relative stability of the injection pressure of Ap rec (p)/p < 
10~ 3 . The performance of the beam tube cooling system and the corresponding 
temperature monitoring has been tested with an almost full scale demonstration setup 
of the WGTS using original components. Fourier analysis of the measured temperature 
fluctuations on the beam tube confirm the functionality of the cooling circuit design 
and in particular prove that temperature or mass flow fluctuations in the primary LHe 
circuit are not transferred to the cooling circuit of the beam tube. 

Complementary to the control and monitoring of the experimental parameters 
defining Af, another method has been developed to monitor fluctuations in M close 
to real-time; it is based on measuring tritium purity and source activity simultaneously. 

The tritium purity, which is also of interest in its own right, is hereby determined 
with a Laser Raman setup (LARA); LARA has also exceeded its specification and has 
reached 0.3% statistical uncertainty in T 2 monitoring of a low mole-fraction tritium 
sample. Based on this result, it is expected that the system reaches 0.1% stat. 
uncertainty in < 60 s for KATRIN relevant gas compositions. The system was calibrated 
with non-tritated gas mixtures. For each isotopologue, specific response functions 
were determined which vary by up to ~10%. Further improvements of the accuracy 
are expected once quantum mechanical transition matrix elements are included whose 
accurate experimental verification has just been completed. 

For the activity monitoring, three different approaches are pursued. Two methods of 
measuring the source activity can be implemented in the Control and Monitoring System 
(CMS) at the rear end of the KATRIN setup: the /3-electron-current through the rear 
wall of the KATRIN experiment can be measured with a Faraday cup, and Beta Induced 
X-ray spectroscopy makes use of the X-rays produced by these /3-electrons to monitor the 
source activity. The third option, the forward beam monitor detector FBMD monitors 
the edge of the flux-tube in forward direction towards the main spectrometer. All 
three options are on the R&D way to demonstrate their suitability to measure with the 
required precision of 10~ 3 under KATRIN conditions. 

Since most of the systems are now operational and perform beyond expectation, we 
are confident that the challenging control and monitoring requirements for the operation 
of the KATRIN WGTS will be met and even exceeded. As a consequence of this 
better-than-expected performance, it might even be possible to further reduce systematic 
uncertainties. 
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